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1. Motivation and objectives 

In premixed turbulent combustion modeling, the most common approach is based 
on averaged transport equations for statistical moments (mean, variance) of thermo- 
chemical variables. Most of the recent developments axe for systems in the flamelet 
regime and consider infinitely fast chemical reactions compared to the turbulent 
mixing (Bray 1980, Peters 1986). Under this assumption, the turbulent flame can 
be represented as a collection of laminar flamelets stretched and convected by the 
turbulent flow field. Thus, several authors have proposed to describe the turbu- 
lent flame by a transport equation for the surface flame area per unit of volume £ 
(Bideaux et al 1994, Candel Sz Poinsot 1990, Pope 1985-a, Trouve & Poinsot 1994) 
or by a transport equation for an arbitrary scalar field (concentration, temperature) 
G(x*,f) = Go which describes the flame sheet as an infinitely thin interface sepa- 
rating fresh and burnt gases (Kerstein et al. 1988, Peters 1992, Williams 1985-b). 
Another possibility is to derive a transport equation for the dissipation of the fluc- 
tuations of the concentration leading, in the limit of high Damkohler number, to a 
balance equation for £ (Mantel & Borghi 1994). 

Another general approach consists of directly calculating the probability den- 
sity function (pdf) of the thermochemical variables through its transport equation 
(Borghi 1988, O’Brien 1985, Pope 1985-b). This method avoids the modeling of 
the chemical reaction rate and can more easily take into account reduced chemical 
schemes or complete chemical schemes as well (Vervisch 1991). However, the lim- 
itation of these models is the high computational cost of the Monte-Carlo method 
used and the sub-models for the molecular mixing terms. 

Recently, Bilger (1993-a,b) has proposed an alternative method based on con- 
ditioned moment closure (CMC) applied to turbulent reacting flows for both non- 
premixed and premixed systems. This technique goes beyond the zone condition- 
ing methods developed and applied to non reacting flows in the case of a two- 
dimensional mixing layer by Libby (1975) and extended to the turbulent mixing 
of scalars by Dopazo (1977). Like pdf approaches, the CMC method is valid in 
the different regimes of combustion and can take into account reduced or complete 
chemical schemes. Moreover, the chemical reaction term can be decomposed using a 
Taylor expansion (Bilger 1993-b) and can be estimated without modeling. Although 
CMC methods in turbulent reacting flows are of recent origin and still under de- 
velopment, applications to real configurations for non-premixed combustion have 
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been already performed using detailed chemistry (Smith et al 1992). Comparisons 
with experimental data obtained for a turbulent jet-flame show encouraging results 
(Smith et al. 1993). As in classical flamelet models (Bray 1980), the dissipation 
rate of a conserved scalar appears in the CMC approach and is a crucial quantity 
which must be modeled. To validate models for the conditional dissipation rate, 
Mell et al. (1994) use results from direct numerical simulation of a reacting mixture 
in homogeneous turbulence whereas Li & Bilger (1993) utilize experimental results 
in a reactive scalar mixing layer. Conditional statistics concerning the dissipation 
of temperature fluctuations in a non-reacting system are also available (Jayesh & 
Warhaft 1992). 

Nevertheless, for premixed systems only theoretical developments exist (Bilger 
1993-b) and certain problems are still unresolved. This is particularly the case for 
the conditional mean velocity and the conditional mean scalar dissipation rate. This 
last term characterizes the molecular mixing (or the rate of dissipation of the scalar 
fluctuations) and must be treated carefully in turbulent premixed combustion as 
has been shown in Borghi (1990) and Mantel & Borghi (1994). 

But, unlike diffusion flames, scalar dissipation is a quantity difficult to measure 
in premixed flames because of the weak spatial resolution of experimental devices 
(Bilger 1993-c). Experimental results for the scalar dissipation have not yet been 
published. Since for premixed systems no data axe available, another possibility 
consists of analyzing results of DNS of turbulent premixed flames. Among the 
available databases, we cite the pioneering work of Trouve & Poinsot (1994) who 
have simulated the interaction between a premixed flame and a three-dimensional 
isotropic decaying turbulent flow field with realistic heat release. The analysis of 
these numerical results will provide new information concerning the conditional 
statistics in turbulent premixed flames and eventually lead to new models for these 
unclosed terms. 

The objective of this paper is to briefly introduce CMC methods for premixed sys- 
tems and to derive the transport equation for the conditional species mass fraction 
conditioned on the progress variable based on the enthalpy. Our statistical analy- 
sis will be based on the 3-D DNS database of Trouve & Poinsot (1994) available 
at the Center for Turbulence Research. The initial conditions and characteristics 
(turbulence, thermo- diffusive properties) as well as the numerical method utilized 
in the DNS of Ttouve &; Poinsot (1994) are presented, and some details concern- 
ing our statistical analysis are also given. From the analysis of DNS results, the 
effects of the position in the flame brush, of the Damkohler and Lewis numbers 
on the conditional mean scalar dissipation, and conditional mean velocity axe pre- 
sented and discussed. Information concerning unconditional turbulent fluxes axe 
also presented. The anomaly found in previous studies (Libby & Bray 1981, Moss 
1980) of counter- gradient diffusion for the turbulent flux of the progress variable is 
investigated. 
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2. Accomplishments 

2.1 Equation for the conditional average of the species mass fraction 

The transport equation for the conditional average of the mass fraction of a 
reactive species has been derived by Bilger (1993-b) from the balance equations for 
the species mass fraction Yi for the species i and for the progress variable c based 
on the enthalpy which will be defined below. Its derivation is similar to that used 
for non-premixed flames (Bilger 1993-a). Klimenko (1990) has a different approach 
which has not been implemented in premixed flames. 

The balance equation for the mass fraction Y{ is: 


P 


dY rr 

-ar +fDk 


dx k 


d 

dxk 




a) 


where T>i is the molecular mass diffusivity for the species i, p the mass density, Uk 
the velocity vector, and Wi the chemical reaction rate. 

For the progress variable c, we use the special definition proposed by Bilger (1993- 
b): 
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Here, h is the total enthalpy, h* represents the sensible enthalpy, and the subscript 
u characterizes the unburnt gases. The sensible enthalpy rise A h* is given by the 
adiabatic burnt value less the unbumt value. Thus, the progress variable c varies 
to 0 in the fresh gases to 1 in the fully burnt gases for adiabatic flames. For non- 
adiabatic flames c can go beyond unity, this being considered useful for studying 
NO formation and CO burnt out. 

The transport equation for h and h 9 are respectively: 
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where a represents the thermal diffusivity of the mixture, p the pressure, $r the 
radiative heat flux, and A h{ the heat of formation for the species i . 

According to its definition given in Eq. (2), the equation for c is directly obtained 
by combining Eqs. (3) and (4): 
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Here, the pressure rise and the radiative heat losses have been neglected. 

The conditional average of the species mass fraction Yi conditioned on the progress 
variable c(x*,t) being at a chosen value £ can be defined: 

Qi«,x k J) = (Yi(x k ,t) | c(x k ,t) = C> = (Yi | C) 


( 6 ) 
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From this definition and from Eqs. (1) and (5), Bilger (1993-a) has obtained the 
balance equation for Q,: 


<' I <>£ + <+ 1 l <>g = I C>(^^ 1 0^ + <* 1 0 
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with 
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Here, Le t represents the Lewis number for the species i, u* and y f * are the fluctu- 
ating parts of respectively Uk and Yi about their conditional means, and Q yC is the 
derivative of Q with regards to c. The probability density function of c defined at 
c = C is given by P c (C)- 

If the Reynolds number of the flow is sufficiently high, the molecular diffusion 
fluxes present in eg are negligible compared to the turbulent flux (Bilger 1993- 
a). The conditional turbulent flux e y has been observed negligible in the case of a 
reacting mixing layer (Bilger 1993-a) and can probably be neglected in premixed 
flames. 

The conditional mean chemical reaction rate can be approximated using a Taylor 
expansion for (u>, j £) (Bilger 1993-b). This technique has been used by Smith et 
ai (1992) for the prediction of radicals and NO in a turbulent diffusion jet flame 
and comparisons with experimental data lead to satisfying agreements. 

Even after these approximations, some terms of Eq. (7) remain unclosed and need 
to be modeled. Our attention will be focused principally on the conditional mean 
velocity (Pjt | C) and on the conditional mean scalar dissipation N$ defined by: 
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These two terms are also crucial quantities in models using pdf methods. They 
appear explicitly in the transport equation for the pdf P( C): 


(p I C)§ + (p I C><% I C>f£ = £(<»« I <)p) + ^((P I <)N ( P) (12) 


where w c is the last term appearing on the right hand side of Eq. (7). 
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2,2 The numerical experiment of Trouve £ Poinsot (1994) 

Since measurements in premixed flames are limited due to spatial resolution prob- 
lems (Bilger 1993-c), DNS offers the possibility of having access to specific quantities 
of statistical interest such as unconditional and conditional moments (mean, vari- 
ance, etc). The numerical results used in the present analysis come from the DNS of 
turbulent premixed flames performed by Ttouve and Poinsot (1994). These authors 
have studied the interaction between an initially planar premixed laminar flame 
with realistic heat release and a three-dimensional decaying isotropic turbulent flow 
field. Their code fully resolves the compressible Navier-Stokes equations using a 
sixth order spatial scheme (Lele 1992) and a third order temporal scheme (Wray 
1992). Due to the heat release and the resulting gas expansion, the boundary con- 
ditions in the direction of propagation of the flame are inflow/outflow non- reflecting 
(Poinsot and Lele 1992), whereas periodicity is used in the two other directions. 

The chemistry is described by a single-step irreversible reaction using a classical 
Arrhenius law. Following the notation of Williams (1985-a), the reaction rate for 
the deficient species (subscript JR) can be expressed: 

(i3) 


where: 


0 =(T - T u )/(T h — T u ) is the reduced temperature 
<p =(7* — T u )/T 6 represents the heat release parameter 

(3 —ipEa/ iZ°T& is the ZeVdovich number 

A =£exp(— /3/V) is the pre-exponential factor 

The subscripts u and b represent respectively the unbumt and burnt gases and E a 
the activation energy of the reaction, R° being the universal gas constant. 

Moreover, in the simulations, the temperature dependence of the transport coef- 
ficients is taken into account by a power law: 



where p represents the dynamic viscosity and b is a constant equal to 0.76. The 
thermal conductivity and molecular diffusivity of the reactants are determined by 
considering constant Prandtl and Schmidt numbers. 

Different simulations have been performed for various Lewis number (0.8, 1.0, 
1.2), the other physical parameters being the same for all the configurations. At 
time t = 0, the turbulent Reynolds number based on the Taylor microscale A is 
Re\ = 50 and the initial Damkohler number (defined by Da = (A/u')/(a u /S£ o )) 
is Da — 0.4. Here, a u represents the thermo-diffusivity in the fresh gases, Si 0 the 
laminar flame, speed and u 1 the rms turbulent velocity. These initial conditions 
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FIGURE 1. Evolution of the DNS in the diagram of premixed turbulent flames - 
comparisons with the operating domain of classical experimental devices and S. I. 
engines. 

correspond to the ratios u f /S l q = 10 and lt/6p = 5.2 (It is the integral length scale 
of the turbulence and 6p the laminar flame thickness defined by 6p = oc u j S l 0 )- 

The flames investigated in the DNS evolve in the thickened wrinkled flames regime 
of the Borghi diagram for turbulent premixed flames (Borghi 1985) as is shown in 
Fig. 1. Also schematically shown in this diagram is the operating domain of classical 
experimental devices and S. I. engines (Abraham 1985). 

In the following analysis, all of the cited turbulence quantities are obtained from 
a simulation of 3-D decaying isotropic turbulence with no combustion. All cases 
have been performed using a 129 3 computational domain. 

2.S Statistical analysis 

In our analysis, since radiative heat losses and pressure rise are neglected and the 
specific heat is constant, the progress variable c defined in Eq. (2) is equivalent to 
c = (T — T u )/(Tb — T u ) and is equal to 0 in the fresh gases and to 1 in the fully 
burnt gases. In the calculation of the scalar dissipation N = a(dc/ dxk)(dcj dx\ t), 
the thermo- diffusivity a is temperature dependent according to the relation (14) 
and by using the constant Prandtl number ( Pr — 0.75). 

In the determination of the conditional statistics, the progress variable space £ 
is divided into 100 sub-domains. In each point of the computational domain, the 
variables of interest are calculated and stored in arrays corresponding to the sub- 
domain of £ at this point. Conditional averages are then obtained by dividing the 
variable contained in the arrays by the number of samples of each sub-domain. 
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Since in the DNS the direction of propagation of the flame front is following the x 
direction, only (U\ | C) is of interest and will be investigated in the present analysis. 

For the unconditional statistics, since we have access to only one realization, the 
problem is assumed to be homogeneous in the y — z directions (x being the direction 
of propagation of the flame front). Then, the computational domain is divided into 
slices perpendicular to the x direction, and statistics are calculated for each plane 
through the turbulent flame front. The averaging of a variable Q varying spatially 
and temporarily is then estimated as it is proposed by Trouve & Poinsot (1994): 

(Q(*,t)) = J Q( x ,y,z,t)dydz (15) 

where L y and L z are the size of the computational box in the y and z directions. 

Through the flame brush c, the density weighted mean of c increases monotically 
with x. Consequently, the results will be presented using the mean progress variable 
c instead of the x coordinate. 


2.4 Results 

2.4-1 The conditional dissipation rate N ^ 

Because it represents the small scale mixing or molecular dissipation of a scalar 
(concentration, temperature), the dissipation rate of the progress variable N is a 
crucial quantity in turbulent combustion (Borghi 1990, Bray 1980, Mantel & Borghi 
1994). In the case of turbulent premixed combustion, the mean reaction rate is di- 
rectly proportional to N as has been demonstrated by Bray (1980) in the limit of 
high Damkohler numbers. Often, this term is modeled utilizing the dissipatioirtime 
scale of the turbulent kinetic energy and a suitable constant, so that N oc c /,2 /r t 
(Bray 1980, Spalding 1971). However, further studies have shown that the pro- 
portionality constant depends on at least two parameters describing the structure 
of the premixed flame such as the Damkohler and turbulence Reynolds numbers 
(Mantel & Borghi 1994, Said & Borghi 1988). In those studies, it appears that the 
ratio of the two dissipation time scales R = Tt/r c (r c represents the dissipation time 
scale of the scalar) is a function of the ratio k 1 / 2 /Sl q , which can also be defined by 
k 1 / 2 /Sl 0 oc (Re t /Da) 1 / 2 ( k represents the turbulent kinetic energy). Said & Borghi 
(1988) proposed an algebraic formula for R whereas Mantel & Borghi (1994) have 
derived a transport equation for N. 

For the conditional dissipation rate, few theoretical developments exist, and ex- 
perimental and numerical data are only available for non-premixed combustion or 
non-reacting flows. Both Mell et al. (1994) in their numerical simulations and Li & 
Bilger (1993) in their experimental study have modeled the conditional dissipation 
of a conserved scalar N v using the unconditional one and a frequency related to the 
large scales of the turbulence: 


N V *N Z = RZ " 2 1 


(16) 
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where R is the time scale ratio assumed to be constant, e represents the dissipation 
of the turbulent kinetic energy, Z is the mixture fraction in the physical space, and 
7} the mixture fraction in the domain of phase. These authors conclude that this 
classical closure gives good predictions for their experimental and numerical results, 
although more detailed modeling of is better (Mell et ai 1994). 

However, this kind of closure apparently valid for non-premixed systems cannot 
be applied for premixed flames because the turbulent mixing is not the only cause 
of scalar variance annihilation. Chemical reaction is also involved and has to be 
considered in the modeling of Nq. This is illustrated in Fig. 2 which shows the 
evolution of R across the flame brush. The ratio R varies significantly through 
the flame brush ranging from 0.3 to 1.5. This is consistent with the values found 
by Beguier et ai (1978) in the case of non- reacting flows. Only a slight effect 
of the Lewis number is observed on R (see Fig. 2-a). As the time goes on, we 
note in Fig. 2-b that R is less than unity in the most part of the turbulent flame 
front (0.1 < c < 0.9). This indicates that the turbulence dissipates more rapidly 
than the fluctuations of c. In such a case, it clearly appears that the well-known 
Eddy Break-Up model will significantly overestimate the mean reaction rate since 
Webu CC c(l - o)/r*. 

Fig. 3 shows the form of iV^ am . = iam T ch ( T ch represents a chemical time 
scale defined by r c ^ = a n j S\ o ) in a planar laminar flame front for different Lewis 
numbers. As expected, N$ lam is insensitive to the Lewis number since the progress 
variable is expressed in terms of temperature instead of the concentration of fuel, 
for instance. Unlike Bilger (1993-b), N $ does not have a symmetric shape, but 
has a maximum value towards the burnt gases. This asymmetry is due to the 
high activation energy for the kinetics, whereas Bilger (1993-b) is for a full kinetic 
mechanism. We can note the quadratic behavior of N^ lam in the preheat zone 
(£ < 0.6) as given by the convection diffusion balance: 



which leads directly to N{ Um = (T/T u ) l ~ b ,£ 2 . 

Fig. 4 shows the evolution of (normalized by the laminar value N $ Um ) for the 
turbulent flame with Le = 1.0 versus £ for different times (t/ro ranges from 1.4 to 
4.5, To being the eddy turn-over time). At the earliest time results are shown after 
one eddy turn-over time, and thus the flame has had the time needed to adapt itself 
to all the spectrum of the turbulence. Since in the DNS positive stretch occurs (on 
average) more often than negative stretch (Trouve & Poinsot 1994), the profiles of 
concentration and temperature are steepened, and we can expect higher values for 
than in the laminar case. At earlier times, this enhanced mixing appears to 
be very strong in the preheat zone but also may be dependent on the Damkohler 
number. Since the Damkohler number is low, the instantaneous flame front does 
not respond immediately to the turbulent flow field, and in the preheat zone, is 
2 or 3 times greater than the laminar value. The reaction zone, however remains 
almost unaffected by the turbulent mixing, and we observe an only slight increase 
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FIGURE 2. Profiles of R = Tt/r c versus c. a) for different Lewis number at 

t/ T Q = 4.5: Le = 0.8: ; Le = 1.0: ; Le — 1.2: . b) for different 

times and Le = 1.0: t/ro ~ 1.4 - Ret = 50 - Da = 0.5: ; t/ro = 2.7 - Re t = 33 

- Da = 0.6: ; f/r 0 = 3.6 - i?e* = 29 - Da = 0.8: ; t/ro = 4.5 - iZe* = 25 - 

Da = 1.0: 

of of approximately 20%. As the time goes on, the Damkohler number increases 
due to the decay of the turbulence. Thus, the turbulent flame comes more into the 
flamelet regime and behaves more and more as in the laminar case. At this 
point, there is no distinction between the preheat zone, and the reaction zone and 
N{/Nt larn is almost constant and close to unity in all of the domain of £. 

The variation of N^/N^ larn across the flame brush is represented in Fig. 5 for 
different values of It appears that N^/N^ lam is almost constant through the 
turbulent flame front and is only weakly dependent on c. Thus, conditional statistics 
for can be extracted from all the computational domain rather than in y — z 
planes. 
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FIGURE 3. Profiles of Nq for laminar flame versus (," given by DNS. L e = 0.8: 
: Le = 1.0: ; Le = 1.2: 



progress variable 

FIGURE 4. Profiles of N^/NQ am versus ( for Le = 1.0. f/r 0 = 1.4: ; 

t/r 0 = 2.7: ; t/r 0 = 3.6: ; t/ro — 4.5: . See Fig. 2 for the 

corresponding Damkohler numbers. 

The influence of the Lewis number on N ^ is also investigated and presented in 
Fig. 6 for 2 different times of the interaction. At time t/ro = 1.4 (corresponding 
to Da = 0.5), we notice a significant effect of the Lewis number in the reaction 
zone (C > 0.7). In the preheat zone, no effect of the Lewis number is observed. 
The turbulent mixing predominates and diffusivities (of species and temperature) 
Eire too slow compared to the chemistry to allow the profile of c to behave like in 
the laminar case. In the reaction zone, the end of the temperature profile seems 
to be affected, illustrating an effect of the Lewis number. The gradients of c are 
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FIGURE 5. Evolution of N^/N^ tam versus c for different values of £ at t/ro = 4.5 

and Le = 1.0. C = 0.0: ; C = 0.25: ; C = 0.5: ; C = 0.75: ; 

C = 1.0: 


steeper for Le = 0.8 and flatter for Le = 1.2 (for Le = 1.2, we can note that N ^ 
decreases significantly below the laminar value). There is no clear explanation for 
such behavior. We can argue that for a positive stretch, the local reaction rate 
increases (decreases) when the Lewis number is smaller (greater) than unity and 
affects the end of the temperature profile increasing (decreasing) the gradients of c. 

For time t/ro = 4.5 (which corresponds to Da = 1.0), N ^ is affected uniformly by 
the Lewis number and no distinct behaviors exist between the preheat and reaction 
zones. The gradients of c are steeper everywhere in the ( space for Le = 0.8 and 
weakly flatter (especially in the reaction zone) for Le = 1.2. This result is quite 
surprising because we can expect that when the Damkohler number increases, the 
turbulent flame behaves more and more in the flamelet regime. Thus locally the 
flame should have the structure of a laminar flame, and consequently the 3 curves of 
Fig. 6-b should merge together. This comment concerns mainly the case Le = 0.8 
for which the gradients of c are significantly higher than the two others cases. Since 
in this case the Lewis number is smaller than unity, we can hypothesize that thermo- 
diffusive instabilities have become non-negligible and provoke additional stretch that 
is not present in the cases Le = 1.0 and Le = 1.2. 

In conclusion, from the Fig. 6, it appears that the ratio N^/N^ lam depends 
strongly on the Damkohler number, particularly in the preheat zone. At this point, 
we cannot make a conclusion on the dependency on Lewis number since we are 
not able to separate the effects of the turbulence and the effects of thermo- diffusive 
instabilities. 

In order to get more information on the role of the mixing on the structure 
of the flame, the product P((^;c).N^ is extracted from the DNS for Le = 1.0 at 
t/ro = 4.5. Fig. 7 shows P((;c).N^ (non-dimensionalized by r c h) as a function of ( 
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0 0.2 0.4 0.6 0.8 1 



Figure 6. Profiles of N^/N^ am versus (. a) t/r 0 - 1.4 - b) t/r 0 = 4.5 Le = 0.8: 
; Le = 1.0: ; Le = 1.2: 


at different locations across the flame front. It appears that the dissipation occurs 
almost everywhere in the progress variable space, but most of it takes place in the 
reaction zone associated with values of £ comprised between 0.6 and 0.9. Higher 
dissipation occurring in the preheat zone is only observed in regions of the turbulent 
flame located near the fresh gases. 

The density- weighted pdf P(£) is shown in Fig. 8 for different sections of the 
flame brush at different instants and Lewis number. For all the cases studied here, 
the pdf P(C) exhibits a strong bimodal shape corresponding to fresh and fully burnt 
gases. However, P( () is not purely bimodal and a broad peak on the fresh gas side 
(which characterizes regions of the flame brush where instantaneous flame fronts 
are thickened) is observed. Between these two peaks (for 0.3 < £ < 0.9), the shape 
of P(C) is rather flat and ranges from 0.2 to 0.8. This represents intermediate 
events occurring in the instantaneous flame fronts. This behavior can be expected 
according to the schematical description given in Fig. 1 concerning the regime of 
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FIGURE 7. Evolution of P(Q.N{ (non-dimensionalized by r c *) versus £ across the 

flame front at t/r 0 = 4,5 and Le = 1.0. c = 0.05: ; c = 0.25: ; c = 0.5: 

; c = 0.75: ; c = 0.95: 

propagations of these flames. 

Since most of the available data for turbulent premixed combustion axe obtained 
using unconditional statistics, it is of interest to calculate unconditional variables 
from their conditional values. The unconditional statistics for a variable Q can be 
directly obtained from the conditional value (Q | C) and the probability density 

function P(£): 

Q= f {Q I C)P(O d C (is) 

JO 

if (Q I C> is known, it appears that the pdf P(C) has to be carefully estimated. 
As we said in the introduction, the pdf P(C) can be directly obtained by solving 
its transport equation (an expression for P(() can be found in Borghi 1988, Pope 
1985-b) or estimated by assuming realistic shapes for P(C) (Borghi 1988, Bray 1980) 
and predicting only its first two moments, c and c" 2 . 

A way to estimate P(£) is to choose a (3 function as proposed by Janicka & 
Kollmann (1978): 

p( 0 


— 1^ and b = Fig. 9 represents the comparison be- 

tween the pdf P(C) given by the DNS and the pdf estimated from Eq. (19) for 
different positions across^the turbulent flame brush. The parameters a and b axe 
calculated using c and c" 2 extracted from the DNS. The agreement between the 
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0 0.2 0.4 0.6 0.8 1 



0 0.2 0.4 0.6 0.8 1 

progress variable 


FIGURE 9. Profiles of the pdf P(C) versus ( across the turbulent flame front at 
Le = 1.0 and t/ro = 4.5. a): DNS ; b): /? function, c = 0.05: — — ; c = 0.25: 
; c = 0.5: ; c — 0.75: ; c = 0.95: 


DNS and the model of Eq. (19) is relatively satisfying for each of the five sections 
studied across the turbulent^flame front. The (3 function constitutes an interest- 
ing approach to estimating P(() in a presumed pdf approach applied to turbulent 
premixed combustion. 

As we have already mentioned, the pdf P(C) appears strongly bimodal, but a 
significant fraction of events occurs between the two spikes located at C = 0 and 
C = 1. This is probably a consequence of the low Damkohler number. Thus, the pdf 
present between C = 0 and £ = 1 can become a crucial quantity in the estimation of 
some terms which are equal to 0 for ( = 0 and ( = 1. This is particularly the case for 
the conditional dissipation rate for which we have | ( = 0) = | £ = l) = 0. 
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The unconditional dissipation rate defined by: 

n= /VcK^CMC (20) 

JO 

is strongly dependent on the value of P(() between the two spikes, especially when 
Da is low as in this DNS. 

In order to estimate how far we are from a purely bimodal combustion regime, 
we have also studied the fluctuation level parameter g defined by Bray (1980): 


9 = 



( 21 ) 


This parameter varies from 0 for no fluctuation to 1, the maximum possible value 
corresponding to a purely bimodal pdf. In this last case, only fully burnt or fresh 
gases are present in the flow. 

The results in Fig. 10 show that even if the level of fluctuations is high, g is 
significantly below the maximum value of 1.0. Thus, in some regions of the flow, 
intermediate states occurs quite frequently as has already been shown in Fig. 10-a 
where the pdf shows a plateau of about 0.5 when c = 0.5. 

Because different models for turbulent premixed combustion relate directly the 
mean reaction rate to c(l — c) and a suitable frequency, this result is important from 
a modeling point of view. In cases such as this, these models will overestimate the 
mean reaction rate. This is particularly the case in the Eddy Break-up formulation 
(Spalding 1971) and for the BML approach (Bray & Libby 1986). An alternative 
possibility is to directly calculate c" 2 from its transport equation. However, this 
approach requires a good model for the dissipation rate N, which is not an easy 
task. We have already discussed the difficulties and the different possibilities of 
estimating this quantity (Mantel & Borghi 1994, Said & Borghi 1988). 

2-4-2 The conditional mean velocity 

Like the conditional dissipation rate the conditional mean velocity (Uk | Q 
appearing in Eq. (7) has to be modeled. Since very little is known concerning 
( Uk | C), information from experimental data or from DNS is of interest in order to 
propose realistic closure for this term. 

The conditional mean velocity (Uk | £) can be obtained from the conditional pdf 
P(Vk\ C) of U k and c: 

(Uk 10= f°° V k P(V k I c) <*C dV k ( 22 ) 

J — oo 

where V* represents the stochastic vector variable related to the velocity vector 
Uk(xk,t). According to Bayes theorem, P(V * I C) is obtained from the joint pdf of 
V k and C P(V k ; C) and the pdf for C P( C): 


P(V k \0 = P(V k ;0/P(0 


(23) 
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FIGURE 10. Variation of the fluctuation level parameter g versus c. Same legend 
as in Fig. 2. 

The joint pdf of V* and £ can be modeled by presuming realistic shapes for P (Borghi 
&; Dutoya 1978) or by calculating its transport equation (see Pope 1981-1985 for 
the derivation and the modeling of this equation). 

The conditional mean velocity can be extracted from the DNS for the different 
available cases. Since in the y and z directions the conditional mean advection 
term is zero from homogeneity, only the conditional mean velocity in the direction 
of propagation of the flame ( x direction) is investigated. 

Fig. 11 represents the variation of the conditional velocity (U\ | £) across the 
turbulent flame front for different values of £. During the simulation, the velocity 
conditioned on the fresh gases is higher than the one conditioned in the burnt 
gases. Thus, the slip velocity (which is the difference between the mean velocity 
conditioned in the burnt gases and the mean velocity conditioned in the fresh gases) 
is always negative and indicates that the flame displacement is strongly correlated 
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mean progress variable 

FIGURE 11. Profiles of (l^i I C) versus c. a) Le = 1.0 - t/ro = 

t / t 0 = 4.5 c) Le = 0.8 - t/r 0 = 4.5. C = 0.0: ; C = 0.25: 

; C = 0.75: ; < = 1.0: 
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with the velocity fluctuations in the unbumt gas. Thus, on most part of the flame 
front, the burnt gases are animated with negative velocities (toward the fresh gases) 
whereas the fresh gases move with positive velocities (toward the burnt gases). In 
the turbulent flame front, fingers of fresh gases appear in the burnt gases (and 
conversely). 

We can also notice only a slight acceleration of the gases through the laminar flame 
front. All the isopleth contours are moved with the same displacement velocity. This 
is probably due to the fact that the pressure difference across the instantaneous 
flame front is insufficient to produce significant acceleration of the fluid through 
the instantaneous flame front in the x direction. Acceleration of the burnt gas 
occurs in the mean pressure gradient across the turbulent flame brush as a whole 
rather than in the instantaneous flame fronts. It seems that the instantaneous flame 
fronts adjust themselves to provide dilatation without acceleration. If this finding is 
sustained in other studies, it could be an important hypothesis on which to base a 
new approach to the theory of turbulent premixed flame propagation. It is evident 
that the variation of {U\ \ () through the flame brush will be significant in CMC 
modeling. 

2.^.3 Turbulent diffusion fluxes 

In turbulent combustion modeling, the turbulent fluxes of the mean progress 
variable axe classically closed using a gradient transport approximation: 

'TT'n _ dfc 
Uk ° Sc t dxk 

where v% is the eddy viscosity and Sc t a Schmidt number. From previous studies, 
it is known that such modeling is not always sound for turbulent premixed flames 
(Borghi & Dutoya 1978, Libby & Bray 1981, Mantel & Borghi 1994, Moss 1980). 
Both Dutoya &: Borghi (1978) and Mantel et ai (1993) found that the Schmidt 
number Sc t is not constant and varies significantly across the turbulent flame front. 
Moreover, the modeling of the turbulent flux strongly affects the structure of the 
turbulent flame brush. In their study, Mantel et ai (1993) show that the tur- 
bulent flame thickness is over-estimated by about 50% by the gradient transport 
approximation. Libby & Bray (1981) have even formulated a model for counter- 
gradient diffusion due to the large variation of density occurring in reactive flows 
in accordance with the experimental results of Moss 1980). With a sufficiently 
high heat release (r = Tb/T u — 1 > 3), the authors show that a pressure gradi- 
ent present in the flow can preferentially accelerate the low density gases and then 
create counter-gradient diffusion processes. However, the pressure gradient is not 
the only candidate able to provoke counter-gradient diffusion. In a recent study, 
Trouve et ai (1994) show that while counter-gradient diffusion is present at low 
turbulence, in the presence of strong turbulence the flame brush tends to become 
more and more wrinkled and thickened, leading to classical gradient diffusion trans- 
port. The authors also show that the wrinkling of the turbulent flame brush is a 
crucial parameter controlling gradient and counter-gradient diffusion. 
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FIGURE 12. Profiles of u" c” versus c. Same legend as in Fig. 2. 

From the DNS, we have computed the turbulent fluxes u" c” and u"c" 2 across 
the flame brush at different times of the interaction and Lewis number (see Figs. 12 
and 13). For the same heat release as in Libby fc Bray (1981), no counter-gradient 

diffusion is observed for u" c" from the DNS. Similarly, no counter-gradient diffusion 
is observed for the turbulent flux u" c" 2 . 

3. Conclusion 

Conditional statistics for the scalar dissipation N £ and the mean velocity (U\ | C) 
have been extracted from direct numerical simulation. The simulations used in this 
study have been performed by Trouve & Poinsot (1994) and describe the interaction 
between a premixed flame with realistic heat release and an isotropic decaying 
turbulent flow field. The turbulence intensity is at the high end of the range of 
interest and the Damkohler numbers at the low end. 
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mean progress variable 


FIGURE 13. Profiles of u"c" versus c. Same legend as in Fig. 2. 

The conditional scalar dissipation behaves differently in different parts of 
the instantaneous flame fronts. In the reaction zone, is weakly affected by the 
turbulent mixing whereas in the preheat zone, the local gradients of the progress 
variable ( are steeper due to the stretch effects which tend to increase N$. This is 
particularly notable when the Damkohler number is low. As the time goes on, Da 
increases (due to the decay of the turbulence in the simulations), and locally, the 
instantaneous flame fronts resemble more the laminar flamelet. 

As expected from the heat release, the conditional mean velocity (x- component) 
(U\ | C) evolves almost linearly through the turbulent flame brush. We also notice 
a negative slip velocity, which characterizes the wrinkling of the turbulent flame. 
We found that in this simulation, burnt gases are moved towards the fresh gases 
and vice versa, producing a highly wrinkled flame as can be observed in Trouve &: 
Poinsot (1994). Due to a weak local pressure gradient, only a slight acceleration of 
the gases across the instantaneous flame fronts is observed. 
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The weighted pdf P(£) derived from the DNS shows a strongly bimodal shape 
(for ( = 0 and £ = 1). Between these two peaks, non- negligible values for P(Q 
axe observed (typically 0.2) which characterize intermediate states. A satisfying 
agreement is observed between P(() estimated from a /? function and P(C) extracted 
from the DNS. . 

The turbulent fluxes u[c and u[c' 2 across the flame brush are derived from the 
DNS. Contrary to previous studies (Libby k Bray 1981, Moss 1980), no counter- 
gradients have been found in these simulations. According to the recent study of 
Trouve et ai (1994), the presence of a strong wrinkling of the turbulent flame brush 
is more favorable to gradient diffusion. 

The present study constitutes a first analysis of conditional statistics in premixed 
turbulent flames. Further experimental and numerical studies have to be done 
in order to provide new information, notably at higher Reynolds and Damkohler 
numbers. 
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